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ABSTRACT 

Detection of the Warm-Hot Intergalactic Medium (WHIM) using Sunyaev-Zeldovich effect (SZE) 
surveys is an intriguing possibility, and one that may allow observers to quantify the amount of 
"missing baryons" in the WHIM phase. We estimate the necessary sensitivity for detecting low 
density WHIM gas with the South Pole Telescope (SPT) and Planck Surveyor for a synthetic 100 
square degree sky survey. This survey is generated from a very large, high dynamic range adaptive 
mesh refinement cosmological simulation performed with the Enzo code. We find that for a modest 
increase in the SPT survey sensitivity (a factor of 2-4), the WHIM gas makes a detectable contribution 
to the integrated sky signal. For a Planck-like satellite, similar detections are possible with a more 
significant increase in sensitivity (a factor of 8-10). We point out that for the WHIM gas, the kinematic 
SZE signal can sometimes dominate the thermal SZE where the thermal SZE decrement is maximal 
(150 GHz), and that using the combination of the two increases the chance of WHIM detection using 
SZE surveys. However, we find no evidence of unique features in the thermal SZE angular power 
spectrum that may aid in its detection. Interestingly, there are differences in the power spectrum of 
the kinematic SZE, which may not allow us to detect the WHIM directly, but could be an important 
contaminant in cosmological analyses of the kSZE-derived velocity field. Corrections derived from 
numerical simulations may be necessary to account for this contamination. 

Subject headings: cosmology: theory-galaxies:clusters:general-cosmology:observations- 
hydrodynamics-methods:numerical-cosmology:cosmic microwave background 



1. INTRODUCTION 

A significant fraction of baryons in the universe 
are expected to be in the Warm-Hot Intergalactic 
Medi um (WHIM) . Numerical cosmological sim ulations 
(e.g., iDave et al.l l2001l : ICen fc Ostrike^ Il999t) predict 
that somewhere around 30-40% of all baryons exist in 
this phase, characterized by electron temperatures be- 
tween 10^ and lO^K. This shock-heated gas primarily is 
located in large-scale filaments and sheets in the cosmic 
web, strung between clusters of galaxies. Observations 
of this gas (if coUisional ionization is the correct inter- 
pretation) ha ve bee n made in Lyg an d OVI absorption 
(|TriDD et all IMiat iDanforth & ShuIH l2005l . l20f)8h . and 
also are likely detectable in U V absorption lines a nd X- 
ray Hues of OVII and OVIII (jNicastro et all [20051) . De- 
tection of these so-called "missing baryons" is critical to 
understanding the matter content of the universe. 

Recently it has been suggested that the WHIM 
gas may be det ectable using the Sunyaev-Z el dovich 
Effect (SZE) (iH ernandez-Monteagud o et aD 120061 : 
iHallman et all 12007; Afshoxdi_et_al., 200"$. The SZE is 
a consequence of cosmic microwave background (CMB) 
photons being inverse Compton scattered to higher 
energy by hot electrons. This upscattering creates a low 
frequency decrement in the CMB, and a correspond- 
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ing higher frequency increment (for review, see e.g . , 
lRephaelilfT995t lBirkinsha ^'1999'. Car lstrom et al.l l200Zl. 
The SZE is most typically discussed as it relates to 
observations of clusters of galaxies, since the very hot 
(~ KfK) gas there is an effective source of high velocity 
scattering electrons. Detections and analysis of known 
galax y clusters have been made by SZE telescopes 
(e.g., iLaRoaue et "al| 120061: iHalverson eFall HOM) and 
a recent blind search using the South Pole Telescope 
SPT) has detected previo usly unknown galaxy clusters 
Staniszewski et aD 120081 ). Investigators have at- 
tempted to estimate t he SZE signature of galaxy groups 
(|Moodlev et al.l [20081 ) and have made weak detections 
of su p ercluster gas in Corona Bq realis (jBattistelli et al.l 
120061 : iGenova-Santos et all l2005l ). In addition, much 
work has been done to study the possibility of using 
the kine matic SZE to d etermine both clust er peculiar 
veloc i ties (iRep haeli fc Lahav 1991: Haehnelt fc TegmarkI 
Il996t lAgha nim et al. 20 Ql[) and the overall bul k velocity 
field fe.g.. [Kashlinskv fc Atrio-Barandelall2000[ ). 

Recent work suggests that filaments in large-scale 
structure, because of their long path lengths and cor- 
related velocities, may have a unique angular p ower sig- 
nature in the SZE (jAtrio-Barandela et al]|2008[ ). In this 
paper, we show that the SZE may contribute a detectable 
signal from the WHIM gas, observable with extensions 
to the sensitivity of current and upcoming surveys (e.g., 
SPT, Planck Surveyor) over relatively small angular ar- 
eas. We examine the contribution of the filamentary 
WHIM to the angular power spectrum of both the kine- 
matic and the thermal SZE. 

1.1. The Relative Strengths of Thermal and Kinematic 

SZE 
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The SZE can be separated into two spectrally distinct 
effects, one which results from the scattering of the CMB 
by the electrons at the thermal velocity in the gas, and 
the other resulting from the bulk line-of-sight velocity of 
the gas (kinematic SZE). Typically for the purposes of 
clusters of galaxies, the kinematic SZE (kSZE) can be 
considered a small perturbation on the full SZE signal. 
However, in the cooler WHIM gas the relative strength of 
the kSZE increases, becoming comparable to the thermal 
SZE (tSZE). The reason for this is the different depen- 
dency of the kSZE and tSZE on the physical properties of 
the gas. Indeed, the effects should be comparable when 
the bulk radial velocity of the gas is roughly equal to 
the average thermal velocity of the electrons in the gas. 
The thermal SZE signal, characterized by the Compton 
y parameter 

y = [ UeTedl (1) 

depends linearly on Ue and Tg (i.e., is proportional to the 
thermal pressure), while the kinematic SZE, character- 
ized by 

(7T f 

b = / VrTLedl (2) 

C J 

depends linearly on electron density and the line-of-sight 
pecuhar (radial) velocity of the gas. For the kSZE, a 
negative radial velocity (toward the observer) results in a 
positive AT/T. For the purposes of this work, we study 
the SZE signals at 150 GHz (2.1 mm), the frequency 
where the thermal SZE signal has its maximum decre- 
ment against the CMB. For the kSZE, the relative tem- 
perature shift, AT/T, is independent of frequency, and 
is equal to the value of b. The tSZE value for AT/T is 
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which at 150 GHz is equal to —y. Here x = hv/kT. 
The thermal SZE is subject to relativistic corrections, 
which can be importa nt for very hot (> IQkeV) clusters 
(|Itoh fc Nozawal[200l . 

If we take the ratio of 5/y, the relative strength of the 
SZE temperature perturbations at 150 GHz, for a given 
parcel of gas we get 



(4) 
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This expression can be written as 
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A similar calculation is shown in lBirkinshawl ()1999D . For 
a radial velocity of 100 km/s, and a gas temperature of 
5 X 10 K (a temperature typical of clusters), this ratio 
is around 4 percent. But in a gas having the same ra- 
dial velocity but a reduced temperature, say a value of 
T = 10^ K, the ratio is roughly 2. Therefore, for million 
degree gas with a bulk velocity of 100 km/s along the 
line of sight, the magnitude of the kSZE component is 
twice as large as the tSZE. One caveat is that kSZE can 
be either positive or negative, thus could cancel out the 
tSZE in some instances. It is possible for the kSZE to be 
zero when the gas has no peculiar linc-of-sight velocity. 
In addition, gas in cosmological filaments is of relatively 




Fig. 1. — Light cone projection for the low-density WHIM (LDW) 
for the thermal SZE. LDW is defined as lO'^ < T < lO'^K and 
< (5 < 50. Field spans 10x10 degrees. 

low density. However, this gas tends to be more spatially 
extended, thus there can be very long path lengths, cre- 
ating a stronger SZE signal. We should note that AT/T 
for the tSZE is frequency dependent, and on the decre- 
ment side at 30 GHz (where interferometers like CARMA 
operate), AT/T is twice the value of y, while the kSZE 
has no frequency dependence. Thus at lower frequencies, 
the ratio of AT/T between kSZE and tSZE is a factor of 
2 lower than the b/y indicated above. 

I t is also perhaps im port ant to note th a t calc ulations 
bv lFox fc Loebl (fl99l and lYoshida eTaH (l200l for ex- 
ample, show that in hot gas, the electron temperature 
may have a non-negligible equilibration time with the 
(hotter) temperature of the ions in the post-shock regions 
in filaments. A systematically low electron temperature 
impacts the thermal SZE signal directly, as it is linear 
with Tg. The kSZE is unaffected, as it is independent of 

In any case, it is intriguing that the kSZE may serve to 
boost the tSZE signal of the WHIM gas significantly, per- 
haps aiding in its detection in a statistical sense through 
all-sky surveys at the SZE decrement maximum. While 
it is sometimes advantageous to separate the kSZE and 
the tSZE by their spectral signatures in order to measure 
each effect individually, here we also attempt to take ad- 
vantage of their (sometimes) additive effect in order to 
aid in the statistical detection of the WHIM gas. 

2. SIMULATIONS 

The main results shown in this work ar e from analysis 
of the Santa Fe Lightcone, described in iHallman et al.l 
()2007( ). Briefly, this N-body/hydro cosm ological simula,- 
tion was carried out with the Enzo (O'S hea et al.|[2005l ) 
code for the express purpose of generating 100 square 
(10 X 10) degree synthetic sky surveys in the SZE and 
X-rays. The calculation is initialized at z=99 using a 
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Fig. 2. — Light cone projection for the LDW for the kinematic 
SZE. Same field as figure [T] 



lEisenstein fc Hul (|1999D power spectrum with n = 1. 
The cosmological parameters are = 0.3, il-t = 0.04, 
QcDM = 0.26, = 0.7, h = 0.7, and ag = 0.9. We 
suTiulated a cubic volume of dimension 512h~^ comoving 
Mpc with a root grid of 512'^ computational zones and 
seven levels of adaptive refinement using the criteria of 
baryon and dark matter density. We use 512^ dark mat- 
ter particles, resulting in a dark matter mass resolution 
of ^7.3x10^° Mq. 

The light cones are made by stacking 27 slices of the 
simulation data, correcting for their angular scale as a 
function of their physical size, and taking random pro- 
jection axis selections and projection translations to get 
multiple realizations of the light cone for the given cos- 
mology. These slices are projections at various redshifts 
spanning the range 3.0 > z > 0. For this work we use 
100 realizations of the light cone with random projection 
axes and shifts. The final light cone images for this work 
have image resolution of 4096^ for an effective resolu- 
tion of 8.8" per pixel. From our previous work, we have 
shown that this simulation effectively resolves the halo 
mass func tion down to r oughly 5x10^^ Mq. For more 
details see iHallman et aH ([2007i) . 

2.1. Low-Density WHIM Projections 

For this particular work, we have used projections of 
the thermal and kinematic SZE. We have made light 
cones from projections of the kSZE and tSZE from 
gas within the WHIM temperature range (10^ — 10^ K), 
and also limit the gas overdensity to the regime from 
0< S <50 (where 6 indicates the overdensity (gas -f 
dark matter) ratio with respect to pc). The limit of 50 
serves to remove gas that may be very near clusters. In 
a ACDM cosmology with the currently favored param- 
eter values, this is more than a factor of 3 below the 
virial overdensity. Therefore the gas we are probing is 
definitely in the filaments and sheets (and even voids). 
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Fig. 3. — Image of the Compton y parameter for the thermal 
SZE using the low-density WHIM (LDW) gas as described in the 
text. This projection is made from a slice of the simulation box 
at z=0.15, with a width subtending 10 degrees at that redshift 
(comoving physical width of 77/i~^Mpc), and with a depth of 
R:i 128/i"^Mpc comoving. Illustrates how this range of density and 
temperature capture s the filame nts. Projection was made using the 
YT analysis toolkit IITurkll200"8l . yt.enzotools.org). 

In our tests of overdensity threshold, we find very weak 
dependence of the result when changing the overdensity 
cutoff by a factor of a few. We hereafter refer to this gas 
as the low-density WHIM (LDW). The LDW hght cone 
projections for the thermal and kinematic SZE are shown 
in Figures [1] and [2l One reason for the safety factor in 
selected overdensity is that in any simulation, we have a 
finite mass resolution, which means that low mass halos 
are unresolved gravitationally. Therefore the overdensity 
achieved in the unresolved halos is likely lower than the 
typical virial value. The main purpose for our choice of 
limits, therefore, is to ensure that we are fully in the fil- 
amentary WHIM regime. Figure [3] shows a projection of 
a low redshift slice made with the gas properties defined 
in this way, illustrating how this regime of temperature 
and density captures the filamentary gas. 

3. DETECTING THE SIGNATURE OF THE WHIM WITH 
PLANCK AND SPT 

We first investigate the contribution of the WHIM gas 
to the SZE signal in a sky survey. The image histograms 
for the LDW SZE projections are shown in Figure U) 
These histograms are from the raw light cone images, 
which have a native angular resolution of 8.8". It is clear 
that the addition of the kSZE creates a strong increase 
in the number of pixels above a given flux decrement. 
The analogous plot for the full SZE projections (for all 
gas in the simulations) is shown in Figure [51 In this case, 
the total SZE is dominated by the tSZE, and the kSZE 
is a small perturbation on that signal. These figures il- 
lustrate the argument in Section 1.1 very nicely, i.e. the 
kSZE is a much more important contributor in the lower 
density, cooler gas than it is in the hotter cluster gas. 

Next, we add the tSZE and kSZE maps, and smooth 
the light cone maps with a Gaussian kernel where the 
full width half maximum is equal to the beam size at 
the thermal SZE decrement maximum for both Planck 
Surveyor and the South Pole Telescope. We can add the 
two effects for these purposes (adding y and b) since the 



Fig. 4. — Histogram of one of the LDW light cone images in 
Compton y and/or b parameter. Red line indicates the thermal 
SZE, blue is for the kinematic SZE, black is the total of the two, 
an indicator of the total value of AT/T at the thermal decrement 
maximum (Ri 2.1mm/150GHz). 
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Fig. 5. — Same as Figure U but for the full light cone im- 
age including clusters. Red is thermal SZE, black total combined 
signature, and blue is kinematic SZE. Note the relative lack of im- 
portance of the kSZE at the high decrement end. This tail is from 
clusters, where the kSZE is relatively unimportant. 

relative AT/T at the tSZE maximum decrement is equal 
to the value of —y, and the relative temperature shift 
due to the kSZE is independent of frequency. Under the 
assumption that our observation would be made at 150 
GHz, this argument is valid. 

Finally, we add a white noise component with the rms 
equal to the expected sensitivity of the two surveys, make 
a histogram of image values, and attempt to fit a Gaus- 
sian to the histogram. If the signal is below the sensitiv- 
ity of the survey, the fit to a Gaussian should be excellent. 
If the fit is poor, then a signal has been detected. Fig- 
ure [5] shows the mean values for the Gaussian fitting 
for 100 realizations of the LDW-only light cone. It is 
clear that for a four- fold increase in sensitivity in SPT's 
survey, the LDW makes a real contribution to the sky 
signal, and at even a factor of 2 increase (to 5/iK rms 
per beam), there might be a marginal detection. Fig- 
ure [7] shows the result of increasing sensitivity (reducing 
noise) in the full (10 x 10 degree) images using the SPT 
beam size. A sample image histogram for SPT is shown 
in Figure [51 

For Planck, non-Gaussianity appears marginally at a 
factor of 8 increase in sensitivity (0.75 fiK per beam). 
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Fig. 6. — The mean value of from a Gaussian fit to the image 
histograms for 100 light cone realizations with error bars indicating 
the range for 90% of the realizations plotted against the factor by 
which sensitivity has been increased (i.e. noise rms reduced) over 
the expected value for each survey. A value greater than 1 
is indicative of a poor fit to a Gaussian, and that the filaments 
contribute signal to the maps. Blue is for Planck, black for SPT. 

and is clearly there at a factor of 10 increase. The differ- 
ence between these two surveys which favors SPT is the 
angular scale of the beam. SPT has about an arcminute 
beam at this frequency, while Planck's is 7.1'. The scale 
of the filaments (of order Mpc), which is similar to the 
cluster scale, is a better match to the SPT beam size for 
most redshifts. A sample image histogram from an im- 
age smoothed to Planck resolution with its Gaussian fit 
overlaid is shown in Figure [H 

We again note that our simulations assume tempera- 
ture equilibrium between the ions and electrons, which 
is probably not achieved on short timescales in this gas. 
This result should reduce the thermal SZE signature by 
some amount which is uncertain. A survey would have 
to be at least as sensitive as what we have described in 
order to have WHIM contribution to the tSZE signal. 

3.1. Methods for Detecting WHIM Gas in SZE Sky 
Surveys 

While it appears clear that at a given survey sensitiv- 
ity the WHIM gas in filaments (LDW) contributes signal 
above the noise, it is not initially obvious how such a mea- 
surement would be made. In the real universe, we do not 
have maps of the filaments only, and the lack of redshift 
diminution in the SZE results in a significant amount of 
confusion in the survey fields. Therefore, separating the 
contribution of WHIM from that of clusters and groups 
of galaxies may seem an insurmountable challenge. 

However, it is possible that filaments have a unique an- 
gular power signature, such that their contribution could 
be identified separate from clusters. There are significant 
difficulties with this approach, as illustrated in Figures 
[TOl and [TlJ These figures show the relative power for a 
full simulated sky map and for the low-density WHIM 
only projections. Figure fTOl shows this difference for the 
tSZE, Figure[TT]for the kSZE. While there is a significant 
amount of gas in the WHIM phase in simulations, they 
contribute very little to the angular power spectrum, par- 
ticularly in the tSZE. The angular power at the peak is 
orders of magnitude lower in the WHIM for the tSZE, 




Fig. 9. — Gaussian fit to the image histogram for WHIM light 
cone with Planck beam and a factor of 8 increase in sensitivity 
over that expected in Planck's large survey. This is equivalent 
to a AT rms of 0.75 /xK per beam. 



Fig. 7. — Four images of the full 10x10 degree field of the LDW- 
only light cones in total AT/T smoothed with a 1' Gaussian to sim- 
ulate SPT's beam. Upper left: Image with lO/iK rms white noise 
Upper rightiimage with 5/jK rms, lower left: image with 2.5/iK 
rms, lower right: 1.25/iK rms. 
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Fig. 8. — Gaussian fit to the image histogram for WHIM light 
cone with SPT beam and a factor of 4 increase in sensitivity 
over that expected in SPT's large survey. This is equivalent to 
a AT rms of 2.5/^K per beam. 



and more than a full decade lower for the kSZE. More 
troublesome is that the shapes of the power spectra are 
very similar. We should note that for the kSZE, some 
small changes in the stacking algor ithm must be made 
(compared to iHallman et al.l (|2007f )) in order to main- 
tain continuity of structures along the projection depth. 
This means that no shifting in light cone stacking takes 
place unless the depth in projection is equal to the depth 
of the full simulation box. 

Indeed, at these large angular scales the SZE is 
swamped by the primordial CMB anisotropics, and the 
WHIM contributes a very tiny fraction to the power at 
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Fig. 10. — Angular power spectra of the thermal SZE from 100 
light cone realizations with ranges (dotted) within which 90% of 
our realizations fall. Black is for projections of the full simulation, 
blue is for the low-density WHIM as described in the text. 

these multipoles. What is important to note in these 
plots is that the prediction of previous investigators of 
a bump in t he angular power spectrum of the kSZE 
at I « 400 (lAtrio-Barandela et all l2008f ) as a result 
of filaments in the low redshift universe is apparently 
absent from our result. This is unfortunate, as this 
bump would create a unique angular signature which 
one could use to ident i fy th e WHIM in filaments. The 
lAtrio-Barandela et al.l (j2008D work uses an analytic for- 
malism to describe the distribution of gas in filaments, 
and the distribution of the filaments themselves. It is 
not clear which assumptions of their model are not re- 
produced in our simulations, and could account for the 
discrepancy. We will explore this further in later sections. 

It is interesting to note, however, that the relative 
power in the kSZE and tSZE as a function of gas phase 
varies dramatically. Figure [T2| shows the power spectra 
for the full light cone maps with tSZE in black and kSZE 
in blue. Figure [Ol shows the same plot for the LDW gas. 
In this gas, the kSZE dominates at the peak and at most 
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Fig. 11. — Same as Figurc[TOl but for the kinematic SZE. Blue 
is the angular power spectrum of the low-density WHIM, black 
is full light cone, dotted are for the 90% range of 100 light cone 
realizations. 
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Fig. 12. — Angular power spectra of the full light cone image 
using all gas in the simulation volume. Black is for the tSZE, 
blue for the kSZE, with 90% variance of 100 light cone realizations 
indicated by dotted lines. 

scales. This result meshes nicely with the description of 
the kSZE/tSZE ratio in Section 1.1. 

3.2. Effect of the Low-Density WHIM on the Power 
Spectrum 

Aside from the obvious question of whether the fila- 
mentary LDW can be directly detected in SZE surveys, 
there is the related key issue of the WHIM filaments 
creating a contaminating signal that might be misinter- 
preted. For instance, shifts in the SZE angular power 
spectrum might result in shifts in fitted cosmological pa- 
rameters (e.g., the value of erg)- Therefore, it is impor- 
tant to understand the WHIM contribution to the power 
spectra. 

To investigate this question, we have taken a ratio of 
our full simulation light cone projection power spectra to 
power spectra of images where the projection of the LDW 
SZE has been removed. In all cases we have smoothed the 
image with a 1' Gaussian to represent the beam of SPT 
at 150 GHz. We have again done this for 100 light cone 




Fig. 13. — Same comparison as Figure [12] but for LDW only 
light cones. Note that the kSZE has more angular power than the 
tSZE for the gas in the LDW phase. 

realizations, and plotted the 90% cosmic variance scatter 
in the following figures. Figure [M] shows this ratio as a 
function of multipole number for the thermal SZE. While 
there is significant cosmic variance at large scales, the 
mean value (solid line) is a change of order 3% maximum. 
At high I, the ratio is of order 1%. For the kSZE, the 
effect is more important. The lines on Figure [15] show 
the mean and 90% scatter of the power spectra ratios. 
The removal of the WHIM lowers the overall power in 
the map by as much as w8% at / w 100 — 400. Though 
declining in contribution to higher I, at I = 1000 the 
effect is still roughly 5%. 

3.3. WHIM Power m the kSZE 

While this effect is not obvious when looking at the 
WHIM and total power spectra of the images sepa- 
rately, the range of I where the LDW contributes to 
the po wer spectrum is broadl y cons istent with the re- 
sults of lAtrio-B arandel a et al.l ()2008[) . Though the cos- 
mic variance using different light cones is quite large, 
with a large sky are a one should see the effect. In 
lAtrio-Barandela et al.l (|2008D . they find a filamentary 
contribution which moves to lower / as z decreases, as is 
expected due to the reduced angular diameter distance. 
Our lowest z defined in these lig ht cones is 0.05 , where 
the fi lamentary contribution in lAtrio-Barandela et al.l 
(|2008D appears to peak between 1500 < I <2000. We 
see the change at lower where they see the contribu- 
tion coming from z < 0.05. We address this issue further 
below. 

We predict a filament ary WHIM angular power sig - 
nature consistent with lAtrio-Barandela et al.l (|2008f) . 
There are some differences in the two types of analyses 
which may be important. First, lAtrio-Barandela et al.l 
(2008) integrate the SZE signature of their filamentary 
models up to relatively high overdensity (300-500), which 
in our simulation is outside the filamentary regime. In 
addition, our light cones calculate the SZE signature 
using a stack of simulation boxes in which the small- 
est redshift we use is 2=0.05. That limit means that 
low redshift filaments, which are largest on the sky, are 
not represented fully by our method. In their model, 
the contribution of filaments is a strong function of red- 
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Fig. 14. — Ratio of power spectra from the light cones made from 
the full simulation box with LDW gas subtracted to those made 
without subtracting the LDW for the thermal SZE. Dotted lines 
indicate the 90% cosmic variance in our 100 light cone realizations. 
While there is large variation at low I, a result of the small volume 
probed at low z, the mean is consistent with very little change due 
to the presence of LDW gas. 
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Fig. 15. — Ratio of power spectra as in Figure [Ml The removal 
of the LDW results in a reduction of the power at low multipoles. 
The effect has a maximum of 8% at / 100 - 400. 

shift, low z being mos t impo rtant. Also, in contrast to 
lAtrio-Barandela et al.l IIQOH) we represent a 10x10 de- 
gree patch of sky, while their prediction is for an all sky 
signature. 

To explore the potential contribution of foreground fil- 
aments, we have added to our light cones a set of 5 images 
spanning the redshifts fr om 0.05 to 0.00 with a Sz of 0. 01 
(simulating the work of lAtrio-Barandela et al.l (|2008f )). 
For this we have used a smaller simulation of a cubic vol- 
ume of dimension 128ft.^^Mpc with 256'^ root grid zones 
and 5 levels of dynamic refinement with refinement crite- 
ria similar to those used in the original simulation. This 
additional simulation is necessary since we have much 
finer Sz spacing, and can then simulate the contribution 
of low redshift filaments to the kSZ images. For these 
low redshift slices this volume is more than sufficient, 
since very little volume is actually sampled in a fixed an- 
gular area of 100 square degrees at z < 0.05. We have 
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Fig. 16. — Ratio of power spectra as in Figure [TSl in black, over- 
laid with same ratio for light cones adding low z slices as described 
in the text (red line). 5 kSZ slices are added from z=0.05 to z=0.0 
with a spacing of Sz = 0.01 to investigate the power contribution 
of near field filaments to the kSZE power spectrum. The blue line 
is the ratio of the two lines, with the dotted line at equality for 
illustration. The low redshift filaments do not substantively alter 
the result. 

generated 5 fully independent realizations of the slices at 
the lowest z range described above, and added them to 
images analyzed in earlier sections. We then repeat the 
analysis of taking a ratio of power spectra with and with- 
out the LDW filaments. The result is shown in Figure 
[TBI and it shows that the addition of low z filaments does 
not substantively change the results. They simply do not 
contribute significant power to the overall sky signature 
in the kSZE. 

4. DISCUSSION 

While it appears clear that the LDW gas in filaments 
contributes a real signal to the thermal and kinematic 
SZE above a certain survey sensitivity, it is unclear how 
one can take advantage of this effect to "detect" the 
WHIM gas. The signal associated with the WHIM in 
the SZE is washed out by the dominant contributions of 
the hotter, denser gas in clusters and groups, and lies 
in a wavenumber regime dominated by primary CMB 
anisotropics. The unique angular signature that the fila- 
ments contribute may make their detection in the power 
spectrum possible. The WHIM filaments also make a 
contribution to the kSZE power spectrum that could re- 
sult in contamination to studies that look to constrain 
cosmology by the kSZE. 

Additionally, it is possible that higher order moments, 
such as the angular bispectrum or trispectrum (spherical 
transform of the three- and four-po int correlation func- 
tions )( e.g., [KogoI&Kom^i^llOOel) of the images, may 
be more promising measurements of filamentary SZE sig- 
natures. If the combined SZE signal of filaments is non- 
Gaussian fe.g. jYoshida et ar.ii200I' ). a three or four-point 
correlation may show a unique signature of filaments. We 
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leave this analysis to later work. 

This simulation has only adiabatic physics and shocks; 
no non-gravitational physics has been included. How- 
ever, in this density and temperature regime, cooling 
times are typically quite long compared to the Hub- 
ble time. Also, what are usually considered to be the 
primary non-gravitational heating mechanisms (stars, 
AGN) in the universe are found preferentially in high 
density regions. Therefore, the adiabatic approxima- 
tion is accurate for this particular problem. However, 
if episodes of galaxy outflow heating have taken place 
in this low density gas (as may be indicated by its non- 
negligible metallicity), then corrections may have to be 
made. Additionally, the spatial resolution in the low den- 
sity gas in this particular simulation is not very high. 
Our root grid simulation zones have spatial extent of 1 
Mpc, and gas at between overdensity of 1 and 50 has 
spatial resolution of between 250 and 500 kpc per zone. 
Also important is our mass resolution, which means that 
low mass (under 5 xlQ^^ Mq) clusters are gravitation- 
ally under-resolved. Therefore, gas which is in filaments 
in our simulation may properly belong in collapsed struc- 
tures like galaxies. It is not immediately clear which di- 
rection these effects would push the results, since higher 
spatial resolution of the filaments would result in higher 
density, increasing the WHIM SZE signal, and higher 
mass resolution should result in less gas in the WHIM 
phase overall. We will study these questions in later pa- 
pers. 

5. SUMMARY 

Low-density WHIM gas in filaments contributes signal 
to 100-1- square degree sky surveys if the sensitivity is 
sufficiently high. This contribution becomes important 
for an SPT-like telescope survey when the sensitivity is 
increased by as little as a factor of 2-4x. In the filaments, 
the kinematic SZE becomes quite important relative to 
the thermal SZE, particularly in the frequency regime of 
maximum thermal SZE decrement (150 GHz) where SZE 
surveys are being performed. In some phases of the gas, 
the kSZE can be the dominant contributor to the SZE 
signal at that wavelength. 
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